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Abstract. The Poincare gauge theory of gravity has a Lorentz connection with both torsion 
and curvature. For this theory two good propagating connection modes, carrying spin-0 + 
and spin-0~, have been found. The possible effects of the spin-0 + mode in cosmology were 
investigated in a previous work by our group; there it was found that the + mode could 
account for the presently accelerating universe. Here, we extend the analysis to also include 
the spin-0~ mode. The resulting cosmological model has three degrees of freedom. We 
present both the Lagrangian and Hamiltonian form of the dynamic equations for this model, 
find the late-time normal modes, and present some numerical evolution cases. In the late time 
asymptotic regime the two dynamic modes decouple, and the acceleration of the Universe 
oscillates due to the spin-0 + mode. 

PACS numbers: 04.50.-h, 04.50.Kd, 98.80.Cq, 98.80.Jk 
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1. Introduction 

This work reports on an extension of a certain cosmological model, based on the Poincare 
gauge theory of gravity (PGT), which was first announced in fl] and then presented in 
considerable detail in 0. In the latter work it was shown that the dynamic Riemann-Cartan 
geometry (with curvature and torsion) could contribute an oscillating aspect to the Universe 
expansion which could account for the present day observed acceleration. Since then two 
new works have appeared analyzing the dynamics of this model and addressing its fit to the 
cosmological observations L2l SJ . These works have already covered many features of the 
original model in considerable detail. Here we wish to first review the results of the application 
of certain theoretical principles to the PGT. That will naturally lead us to a more appropriate 
description and our extension of the original model. 

One of the outstanding successes of theoretical physics in the latter part of the last 
century which led to a much deepened understanding was the recognition that all the known 
fundamental physical interactions, the strong, weak, and electromagnetic — not excepting 
gravity — can be well described in terms of a single unifying principle: that of local gauge 
theory. Although there are other possible gauge approaches, for gravity it seems highly 
appropriate to regard it a gauge theory for the local symmetry group of Minkowski space time: 
the Poincare group (HO. Such a consideration led to the development of the Poincare Gauge 
Theory of gravity (PGT) HE 0013 QUE). The PGT has a priori independent local rotation 
and translation gauge vector potentials: the Lorentz (i.e., metric compatible) connection and 
the orthonormal co-frame; their associated field strengths are the curvature and torsion. The 
space-time then has generically a Riemann-Cartan geometry. Because of its gauge structure 
and geometric properties the PGT has been regarded as an attractive alternative to general 
relativity. 

The theory includes as exceptional cases Einstein's general relativity (GR) with 
vanishing torsion, the Einstein-Cartan theory with non-dynamic torsion algebraically coupled 
to the intrinsic spin of the source, as well as the teleparallel theories — wherein curvature 
vanishes but torsion does not. The generic PGT has, in addition to the metric familiar from 
GR, a connection with some independent dynamics, manifested in both the torsion tensor and 
additional non- vanishing post-Riemannian curvature components. 

Investigations (especially (HI [T3J) of the linearized theory have identified six possible 
dynamic connection modes, carrying certain spins and parity: 2 ± , l^, ± . It is not possible for 
all of the modes to have good dynamics. The possible combinations of well behaved (carrying 
positive energy at speed < c) propagating modes in the linear PGT theory were identified. The 
Hamiltonian analysis revealed the related constraints |[T4|. Then detailed investigations of the 
Hamiltonian and propagation lfT6l [T5l [T71 [T8l concluded that effects due to nonlinearities in 
the constraints could be expected to render all of these cases physically unacceptable except 
for the two "scalar modes", carrying spin-0 + and spin-CT. 

One mode (referred to as the "pseudoscalar" because of its 0~ spin content) is reflected in 
the axial vector torsion. Axial torsion is naturally driven by the intrinsic spin of fundamental 
fermions; in turn it naturally interacts with such sources. Thus for this mode one has some 
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observational constraints [fT9ll20l . Note that except in the early universe one does not expect 
large spin densities. Consequently it is generally thought that axial torsion must be small and 
have small effects at the present time. The other good mode, + , the so-called "scalar" mode, 
is reflected in the vector torsion. There is no known fundamental source which directly excites 
this mode. Conversely this part of the connection does not interact in any direct obvious 
fashion with any familiar type of matter [|2T1l . Hence we do not have much in the way of 
constraints as to its magnitude. We could imagine it as having significant magnitude and yet 
not being dramatically noticed — except indirectly through the non-linear equations. 

Thus the theoretical PGT analysis led to just two dynamic Lorentz connection modes. An 
obvious place where one might see some physical evidence for these modes is in cosmological 
models. The cosmological homogeneous and isotropic assumptions greatly restrict the 
possible types of non-vanishing fields. Curiously, for the connection and torsion there are 
only two possibilities, which reflect precisely the two spin-0 connection modes. The scalar 
+ which gives rise to a special vector torsion which has only a time component, and the 
pseudoscalar, CT mode, which gives rise to an axial torsion which is the dual of a vector with 
only a time component. Hence the homogeneous and isotropic cosmologies are naturally very 
suitable for the exploration of the physics of the dynamic PGT "scalar modes". 

Thus cosmological models offer a situation where a dynamic Lorentz connection may 
lead to observable effects. Here we will not focus on the early universe, where one could 
surely expect large effects (although their signature would have to be disentangled from other 
large effects), and instead inquire whether one could see any effects of the PGT dynamic 
connection in the present day universe. In particular we will here consider accounting for the 
outstanding present day mystery: the accelerated universe, in terms of an alternate gravity 
theory with an additional natural dynamic geometric quantity: a Lorentz connection 0312. 

The observed accelerating expansion of the Universe suggested the existence of a kind 
of dark energy with a negative pressure. The idea of a dark energy is one of the greatest 
challenges for our current understanding of fundamental physics [[221 l23l |24| . Among a 
number of possibilities to describe this dark energy component, the simplest may well be 
by means of a cosmological constant A. Another popular idea is the quintessence field — 
some unusual type of minimally coupled scalar field — which has received much attention 
over the past few years and a considerable effort has been made in understanding the role of 
quintessence fields on the dynamics of the Universe (see, e.g., [|25ll26ll27l ). 

An alternative is to consider some other gravity theory. Which brings us to our specific 
topic: the possibility of explaining the accelerating universe using a well tested alternative 
gravity theory, one well motivated by both geometrical and physical gauge theory principles. 
We explore the possibility that the dynamic PGT Lorentz connection modes can drive the 
acceleration of the universe. As noted above, there are two spin-0 modes which could have 
dynamical behavior. In [f2l it was shown that the the spin + mode can make the expansion 
rate oscillate, naturally having an accelerating expansion in some periods and a decelerating 
expansion at other times. For suitable choices of parameters and initial data the model can 
account for the supernova observations. Here we show that including the CT mode allows for 
an improved matching. 
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Over the years there have been many studies of PGT cosmology, especially by Minkevich 
and coworkers (see, e.g., Il28l |29l 1301 13T1 1321 ). Using various models they found that it was 
possible for the PGT to avoid singularities, account for inflation, and produce the acceleration 
of the universe (as discussed later in Section 3, their mechanism is different from that of our 
dynamic + mode). A comprehensive early survey of the PGT cosmological models was 
presented quite some time ago by Goenner and Miiller-Hoissen 11331 . Although that work 
only solved in detail a few particular cases, it developed the equations for all the PGT cases — 
including those for the particular model we consider here. However that work was done prior 
to the discovery of the accelerating universe, and torsion was thus imagined as playing a big 
role only at high densities in the early universe. More recently investigators have begun to 
consider various models with torsion as a possible cause of the accelerating universe (see, 

e.g., m El EH EH). 

We have taken another step in the exploration of the possible evolution of the Universe 
with dynamic Lorentz connection spin-0 modes of the PGT. The main motivation is two- 
fold: (1) to have a better understanding of the PGT, in particular the possible physics of the 
dynamic spin-0 modes; (2) to consider the prospects of accounting for the outstanding present 
day mystery — the accelerating universe — in terms of an alternative gravity theory, more 
particularly in terms of the PGT. With the usual assumptions of isotropy and homogeneity 
in cosmology, we find that, under the model, the Universe will oscillate with generic choices 
of the parameters. The + dynamic mode in the model plays the role of the imperceptible 
"dark energy". With a certain range of parameter choices, it can account for the current 
status of the Universe, i.e., an accelerating expanding universe with a value of the Hubble 
constant which is approximately the present one. These promising results should encourage 
further investigations of this model, along with a detailed comparison of its predictions with 
the observational data. 

The remainder of this work is organized as follows: We summarize the formulation of the 
PGT in general and our model with scalar and pseudoscalar modes in Sec. 2, and then consider 
the PGT scalar mode cosmological model in Section 3. In Section 4 an effective Lagrangian 
and Hamiltonian for our cosmological model is presented. This is followed by a late-time 
asymptotic expansion in Section 4 in which certain normal modes are identified. Section 6 
includes the results of our numerical demonstrations for various choices of the parameters and 
the initial data along with a comparison with the supernova observations. The implications of 
our findings are discussed in Section 7 and Sec. 8 is a conclusion. 

Throughout the paper our conventions are as follows: The spacetime signature is 
(—,+,+,+) and c = h = 1. The Greek indices, a,[3,j..., are 4d orthonormal (an- 
holonomic) indices, whereas the Latin indices i,j,k... are 4d coordinate (holonomic) 
indices; they both range over 0, 1, 2, 3. On the other hand the Latin indices a, b, c, d are 3- 
dimensional, with range 1,2,3. 
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2. The Poincare Gauge Theory 

In the Poincare gauge theory of gravity (PGT) |7l[8l, the two sets of local gauge potentials 
are, for "translations", the orthonormal co-frame $ a = e a idx l , where the metric is g = 
— #° cg> $° + 5 a b^ a ® $ b , and, for "rotations", the metric-compatible (Lorentz Lie-algebra 
valued) connection 1 -forms T af3 = T^ al3 \dx l . The associated field strengths are the torsion 
and curvature 2-forms 

T a -df + r^A^^V^/Af, (1) 

R af) ._ dr a/3 + T a^ A ^ = _ R *f3 ^ A ^ (2) 

which satisfy the respective Bianchi identities: 

DT a = R a p A DR a p = 0. (3) 

The PGT Lagrangian density is taken to have the standard quadratic form; qualitatively, 

3 (") 6 (") 

Sf[&,T] ~ A-a R + J2 a n T2 + J2 b " R2 > (4) 

n=l n=l 

(n) (n) 

where T and R are the algebraically irreducible parts of the torsion and curvature and A is the 
cosmological constant. The gravitational field equations obtained from varying with respect 
to the respective gauge potentials $ a i, T af3 j have the qualitative form 

■ A (") (")„ A (")„ 

A — a Ga + 2_. a n[DT + T ) + y. ~ source energy-momentum density, (5) 

n=l n=l 

a T + > a n T + ^ b n DR ~ source spin density. (6) 

n=l n=l 

These are, respectively, second order equations^ for fi a i and T a/3 j. In conjunction with the 
Bianchi identities ©, these two equations yield, respectively, the conservation of source 
energy-momentum and angular momentum statements. 

Here, generalizing an earlier work [2], we wish to examine the dynamics of the special 
case describing the two good PGT dynamic scalar modes in a cosmological model. For this 
two spin-0 modes case we should take b n = except for 6 6 — > b + and 6 3 — > —b~ (see 
Appendix A). For more convenient signs we also make the replacement a n — > — A n . 

The gravitational Lagrangian of our model has the specific form 

3 



1 (n) h+ h~ 

-2A + A R - - Y A n T 2 + —R 2 + —E 2 
2^ 12 12 

n=l 



(7) 



| It should be noted that the PGT is obtained from an action containing quadratic curvature terms, but with the 
connection as a variable they yield 2nd order equations. Higher order equations would result from such an action 
if one used the Christoffel connection — or decomposed the Lorentz connection into a Christoffel part plus torsion 
terms and then treated the torsion and the metric as the fundamental dynamical fields. Such a decomposition is 
alien to gauge principles. 
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where k = 8irG, R is the scalar curvature and E is the pseudoscalar curvature (specifically 
E/6 = -R[oi23] is me magnitude of the one component of the totally antisymmetric curvature). 
The cosmological constant has been included both for generality and a comparison with other 
models. 

In detail, the first field equation, obtained from variation with respect to the orthonormal 
frame, has the components 

^9p,v + AqG^ + A\ ' aT v ^ a + -TvapT^ 013 — T Ufia T a + —g^T^Tafa — T a/ 3^T al3 ^j 
+ ^y^ 1 (g^ a T a - V U T, - \g, u T a T^ + A ~ g A (6V a P^V„ - AP a T^\ 
+ 3P a T a ^\, - AT a Pv P^eJ llx - P a P a g, u + 2P fl P u ) - h —{E 2 g^ - 2ER af3lu e a/3 \) 



+ —{AR^ U - R 9fMU )R = (8) 



b+ 
24 1 

where G^ v = R^ v — \g^R is the Einstein tensor and ^ u is the (in general nonsymmetric) 
material energy-momentum density tensor. 

The components of the second field equation, obtained from the variation with respect to 
the connection, can be decomposed into three algebraically irreducible parts: 

(6(A) - A x ) + b+R)T a ^ - b —ET% u e^ Pl = 0, (9) 

b + V II R-^(6m + + b + R)T fl + h-EP IM =0, (10) 

rV^ - ^(6m~ + b + R)P^ - \b~ET,, = kS^ = 0, (1 1) 

(2) (3) 

where = T% = T%, P M = \e^ af3 T v a p = ±e Mi , Q/3 T% are the torsion trace and axial 
vectors, m + = A Q + A 2 /2 and m~ = A + 2A 3 are the masses of the respective linearized 
modes, and is the axial vector spin density which, for simplicity, we have assumed to 
vanish (this should be a good approximation except at high densities such as those expected 
in the very early universe). The 0~ part couples, as indicated here, to the axial spin vector of 

spin- 1/2 fermions, but the + mode does not couple to any known source. 

(i) 

From © we find the general solution T a : p 1 = 0. From (flOl) and CCD) we learn that the 
torsion trace and axial vectors are controlled by the gradients of two functions. This reflects 
their respective spin + , spin 0~ fundamental nature. However, in view of the non-linearities 
of these relations as well as the geometric significance of the respective "potential functions" 
(i.e., they are the scalar and pseudoscalar curvatures), one can see that it is neither possible 
nor appropriate to resolve them to find new, simpler dynamics for two "scalar potentials". 
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3. The PGT scalar mode cosmological model 

For a homogeneous, isotropic FLRW (Friedmann-Lemaitre-Robinson- Walker) cosmological 
model the isotropic orthonormal coframe has the form 

= dt, d a = a(t) (1 + ^kr 2 )- 1 dx a , (12) 

where k = — 1,0, +1 is the sign of the Riemannian spatial curvature. Here we will consider 
for simplicity only the flat k = case (as far as the observations can tell, this appears to well 
describe our physical universe). 

Because of isotropy, for this k = case the only non-vanishing connection one-form 
coefficients are of the form 

r a = * (t) dx\ Y\ = X(t)e a bc dx c } (13) 

where e abc := e 0abc is the usual 3 dimensional Levi-Civita anti-symmetric tensor. From the 
definition of the curvature (0, one can now find all the nonvanishing curvature 2-forms: 

R 0a = $ dt A dx a - XVe a bc dx b A dx c , R ab = Xe ab c dt A dx c + (^ 2 - X 2 )dx a A da; fe .(14) 

Consequently, the scalar and pseudoscalar curvatures are, respectively, 

R = 6[a _1 * + a~ 2 (^ 2 -X 2 )], (15) 
E = Q[a- 1 X + 2a' 2 X^}. (16) 

Because of isotropy, the only nonvanishing torsion tensor components are of the form 

T\o = f(t)5l T\ c = -2 X {t)e\ c . (17) 

From the definition of the torsion CQ) one can find the relation between the torsion components 
and the gauge variables: 

/ = a - 1 (^-d), x = a - 1 X. (18) 

(Note: the variable $ = —3/ was used in the earlier work O.) 

From the isotropic assumption, the material energy momentum tensor must have the 
perfect fluid form. In this work we focus on the late time behavior. Accordingly we assume 
that the fluid pressure can be neglected, so that the gravitating material behaves like dust with 
a density satisfying pa 3 = constant. Also, we remind the reader that, although we expect the 
spin density to play an important role in the early universe, it is reasonable to assume that the 
material spin density is negligible at late times. 

Due to isotropy, the first field equation ((8]) has only two nontrivial distinct components. 
Expressed in terms of the tensorial quantities and the Hubble parameter H = a/a they are the 
"00" piece 

OA 

+ A + ^H 2 - 3m + {H + f) 2 + 3m- X 2 

+ %E 2 ~ b-E(H + f) X + b -R 2 - h -R[{H + f) 2 - x 2 } = -K P , (19) 
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which contains only first time derivatives of the potentials (and is hence an initial value 
constraint), and the "space-space" piece 

- A + ^-[R - 3(H + ff] + (2m+ - m -) X 2 - ^-{2H + 3H 2 ) 

+ \^E 2 - b -E(H + f) X + ^R 2 - + f) 2 - x 2 } = -K P = 0, (20) 

a dynamical equation for a. 

From the components of the second field equation (11011 II) . we obtain 

b + R = 2(b + R + 6m + )f + 2b~E X , (21) 
b~E = 2b- Ef - 2(b + R + 6m~)x, (22) 

which — along with the definition of the curvature scalars (fl~5l) . (fT6l) — are second order 
dynamical equations for the connection coefficients. 

Before we present our more detailed discussion of these dynamical equations equations 
it should be noted that there are some special "non-dynamic effective cosmological constant" 
cases — that is special cases having one constant magnitude field when certain coefficients 
a/o other field components vanish. Such field components would contribute to the dynamical 
equations certain constant terms that would act like effective cosmological constants. Since 
here we are interested in the generic case with both modes dynamic, we only mention these 
special cases briefly and cite other works where they have been considered in more detail. In 
particular [0 discussed the case with vanishing b~, x an d R = —6m + /b + , while [32J has 
considered the case with vanishing b~ , /, and R = —6m~/b + . It should also be mentioned 
that the earlier investigators have tended to decompose the connection into its Christoffel part 
plus some torsion. Using such a decomposition leads to higher order equations — unless one 
takes A 2 = 0, as these investigators were prone to do (see, e.g., [|28l |29l [30l [32l |33). Here 
we do not make such a decomposition of the connection and need not make such a parameter 
restriction — which would in fact render the + mode non-dynamic. 

Instead of considering the three second-order differential equations (12012 1I22[ ). we can 
transform these second-order equations into six first-order differential equations by including 
the definition of H and (fT5l)-(fT6l) for the six unknown variables: a, H, f, x, R, and E 
(where all except a and H are gauge covariant tensor fields). For some purposes this is 
more convenient; in particular first-order differential equations are more suitable for numerical 
calculations. Combining Eqs. (TT9b and (1201) gives 

- 3A 2 (H + 2H 2 ) + m + R - 4A + 6(m + - mr )x 2 = n{p - 3p) = up. (23) 



Eq. (1231) is considered as a first-order differential equation for H. Moreover, Eqs. (1211) and 
(1221) are already first-order differential equations for R and E. After replacing H in Eq. (PT5T) 
by using (1231) . it is clear that Eqs. (fT5l) and (TT6l) give the first-order equations for / and X- With 
a straightforward re-organization, the six equations are 

d =aH, (24) 

H = ^R+ 2(m+ " m ' } x 2 - 2H 2 - ^ - (25) 
3A 2 A 2 A 3A 2 3A 2 
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; ~^' / - 3Hf+ 3A 2 + 3A 2 > (26) 

X = |-(3iJ + 2/) X , (27) 

* = (29) 
with the constraint equation 

A + ^2# 2 + 3m- X 2 - 3m+(/ + H) 2 + ^i? 2 + h -E 2 

- b~E(f + H)x - h —R[(f + H) 2 - x 2 ] = -k P . (30) 

The constraint equation can be used to replace Kp in Eqs. (1251261 ) to give alternative versions 
of these two equations: 

1 4 m + 3 tt2 2m + -m~ 2 m + 2 

H = " T 2 A + ~ 2 H + — ~ X (/ + H) 

+ i {£ fl2 " T fl [ (/ + ff)2 " ^ + £ B ' 2 " rW + M ' <31) 

- wM r2 - [ (/ + Hf - * 2 i + y 2 - rB</ + m • (M) 

These two alternative equations along with the four other first order equations make a closed 
system for the geometric variables which is more practical for numerical evolution. This 
alternative system will be obtained in another way in the next section. 

For a comparison with GR models, the "00" constraint (fT9l can be considered as a 
generalized Friedman equation: 

3A H 2 = n(p + p r ) + A, (33) 

where the effective energy due to the dynamic connection is 

K Pr = 3m- X 2 -Qm + Hf-3m + f 2 + ^R 2 + b -E 2 -b~Ex(H+f)-^R[(H+f) 2 - X %(34) 



Moreover, the "space-space" Eq. (1201) may be considered as a force-balance equation: 
ra\ _.(p + 3 P ) _ K(pr + 3 Pr ) A 
° \aj 6 6 3' V ' 

where the effective pressure due to the dynamic connection is 

kpy = 2m + f + m + f(AH + f) — vnT\ 2 

+ y 2 r2 + y 2 E2 ~ y e ^ h + f) ~ t r[{h + /)2 " x% (36) 

It can here be seen how p r + 3pr < 0, which is indeed possible, could produce an 
accelerated universe. Although these relations are useful for comparison with other models, an 
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examination of the dynamical terms implicit in p r and pr shows that their PGT dynamic nature 
differs fundamentally from that of their GR counterparts. In particular the ratio tu r := Pr/ Pr 
can take on any value and should not be given the usual equation-of-state interpretation. 

4. Effective Lagrangian and Hamiltonian 

Our cosmological model system of ODEs resemble those of a particle with three degrees of 
freedom. One may suspect that they can be obtained directly from a variational principle. To 
achieve such a goal it is natural to consider imposing the homogeneous-isotropic symmetry 
into the field theory Lagrangian density. We note that it has long been known that imposing 
symmetries and variations do not commute in general. However, for GR they are known to 
commute for all Bianchi class A cosmologies 11371 . We conjecture that this is also true for 
the PGT. In particular our k = model is an isotropic Bianchi I (class A) model, so there 
is good reason to be hopeful. Here we show that (at least for our dust fluid model) imposing 
the FLRW symmetry on the PGT Lagrangian does indeed lead us to the same expressions as 
were found from imposing the symmetry on the field equations. 

Imposing the FLRW symmetry on the Lagrangian density © leads to the effective 
Lagrangian 



It should be noted that (for least action) the coefficients of the quadratic kinetic terms (f 2 , R 2 , 
E 2 ) which contain the time derivatives — see the specific expressions for the curvature and 
torsion Eqs. (|15I16I18I) — must be non-negative. 

Now we use this effective Lagrangian (along with the just mentioned expressions for the 
curvature and torsion) to obtain a conserved energy and three second order equations for the 
gauge parameters: the connection coefficients ij), X and the frame/metric scale factor a . We 
also note that the second order Lagrange equations can be rearranged and combined with the 
formulas for curvature and torsion to give exactly the six first order equations obtained from 
the general 4D covariant field equations specialized to the k = FLRW geometry. Following 
this we will then find the associated Hamiltonian equations. 

Making use of the formulas for the torsion and curvature components in terms of the 
gauge variables, (I15I16I18I) . the conserved energy function associated with L c g is found to be 



where q^ = X, a}. This is a combination that we recognize from our earlier analysis, it is 
just the "00 constraint" (1301) ; its constant magnitude is the physical combination —a 3 p, which 
is indeed a constant because of the dust fluid energy-momentum conservation relation. 

Making use of the formulas for the torsion and curvature components in terms 
of the gauge variables (1 1 51 1 61 1 81) we now obtain the Euler-Lagrange equations 



£ La + + l A2 f - 6A 3X 2 + -R 2 + -e 2 

k 2 2 J A 24 24 



(37) 



8:=q 




- A + ^H 2 - 3m+(/ + H) 2 + 3m- X 2 
k [ 2 

+ %E 2 ~ b-E(H + f) X + b -R 2 - h -R[{H + f) 2 - X 2 ] } , (38) 
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d dL cS _ dL cS _ 

dt dqi dqi 

For the ^ equation: 

<1 dL eS _ d ( sA ^ 2 + b^ a 2 R \ = dL^_ = + h+am + b _ aEX + 3A2fl2/ _ (39) 



dt 84/ dt\ 2 J 

This is a second order equation for \P. (Here and below we dropped for simplicity the overall 
factor of k~ x .) It can alternately be rearranged using (fT8l) to give (1281 ), the first order equation 
fori?. 

For the X equation: 

d dL cS d f \- a 2 E \ = dL*s = b - aE ^ _ b+aXR _ 6Aq(iX _ 12AsaX (4Q) 



dt OX dt \2 J dX 
This second order equation for X can be rearranged using (fT8l) into (|29l) . the first order 
equation for E. 

For the a equation: 

d dL eS d 2 dL cS 2 2 2 6+ 2 2 6+ 2 2 



(-3A 2 a 2 /) = =P = A [a 2 i? - 3(* 2 - X 2 )] + -a 2 R 2 - "-R(V 2 - X 2 



dt da dV J ' da 1 v n 24 2 

+ ^a 2 E 2 - b-X^E + -A 2 a 2 f 2 - 6A 3 X 2 - 3Aa 2 . (41) 

This is a second order equation for a. It can be rearranged into a first order equation for /; 
the result is exactly d32l) . the aforementioned alternative to (T26l) obtained by using d30l) . Now 
using (fl~8l) one can calculate 

if = ^(a- 1 ^) - / = -H(H + f) + a" 1 * - /, (42) 
dt 

then using the just mentioned expression for / and ( TT3T ) one gets (f3Tb . Moreover, from ( fT6l ) 
using (fl~8l) one gets the x equation (|27|) . 

It is remarkable that here in these Lagrange equations and the associated conserved 
energy we get (at least for this dust case) exactly the correct equations for our model — without 
including any explicit source coupling! 

We have cast our system into six first order equations for (3D) tensorial quantities, 
equations which are suitable for numeric evolution and comparison with observations. 
However these equations are probably not in the most suitable form for the most penetrating 
analytic analysis. So we here also present the Hamiltonian equations for our PGT cosmology. 

From the above one can introduce the canonical conjugate momentum variables: 

P Q ^ = -^V /; (43) 
oa k 

dL eS 2 {3A , V 



P x = —4^ = —a 2 E. (45) 
8X 2k 

Now one can construct the effective Hamiltonian: 

V cS = Pad + P^ + P X X-L cS = ^(^ + ^ + §)-(^f + * 2 -X 2 )^ 

6a \ An b~ 6+ / V b + J a 
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+ ^P a - -X VP X + — al 2 + 

a k 



—r 1 + A 

2b+ 



-, (46) 



and obtain the six Hamilton equations: 



dH, 



cff 



Pa 

Px 



3Pa 

0H cS 

OP* 
dH cS 

dPx 
dH cS 

da 

dHes 

dH cS 
dX 



3A 2 a 

K 



# 2 X 2 A 

■>/ P * + IT a ' 

3b + a a a b + 



K 



3b a 



-P 



X 



-XV, 
a 



^P a 12A 



3 X 2 



2A l A 



-P a + ~(VP* + P x X), 
a 

l ^l aX + -(P x m -P 9 X). 
k a 



(47) 
(48) 
(49) 

3a 2 2A 

— + -rr-P^(.50) 
k b + 

(51) 
(52) 



This canonical reformulation should be of considerable interest for further studies of this 
model, since the Hamiltonian formulation is the framework for the most powerful known 
approaches for analytically studying the dynamics of a system, including such techniques as 
the Hamilton- Jacobi method and phase space portraits. 



5. Asymptotic Expansion 



At late times in an expanding universe as the scale factor a becomes larger the field amplitudes 
should be decreasing. We can then expect the quadratic terms in (|38l) to be dominant; hence, 
when the cosmological constant vanishes the late time asymptotic behavior of H, f, x, R, and 
E should have a a~ 3 / 2 fall off. So we reparametrize them according to 



H = ha~ 3/ \ f 



ya 3/2 , x 



xa 



" 3 / 2 , R 



ra 



-3/2 E 



ea 



-3/2 



The current Universe corresponds to a 3 / 2 ^> 1. Substituting (1531) into (I24l)-d30l) gives 

a = a~ 



h 



x 



l ' 2 h, 


-3/2 


"2(m+ - 




A 2 


-3/2 


"4^3 2 


3/2 


3 

2yr + - 


-3/2 


—2xy — 


3/2 


3 

2ey+ 2 



-x - 
3 

— / 

2 

2b~ 



Y 



K 



3A 



Poa 3 (0) 



-p a 3 (0) 



-ex 



3A_ 
12m+ 



rrv 
_l r 

3A 2 
A 



3A, 



r + a 



3/2 



4A 
3A 2 
4A 



3Ao 



-hx 



+ 



26 H 



-rx 



12m" 



-x. 



- a 3 np = -Kp a 3 (0) = a 3 A + 



^/i 2 + 3m"x 2 + b —e 2 + U —r 2 - 3/// + (// - // V 



6+ 
24' 



(53) 

(54) 
(55) 

(56) 

(57) 

(58) 

(59) 
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b + 

b~ex(h + y) + —r((h + y) 2 — x 2 ) 



(60) 



Now let us restrict our further considerations to the vanishing cosmological constant case. 
With a 3 / 2 1, Eqs. (I54l)-(l59l) show that the spin-0 + mode and the spin-CT mode are 
becoming decoupled and the asymptotic Hubble rate h is influenced purely by the spin-0 + 
mode. Dropping the higher order terms (with A = 0), gives the energy constraint 

- a 3 Kp = -Kp a 3 (0) = —h 2 + 3m-x 2 + — e 2 + —r 2 - 3m + (y + h) 2 (61) 

2 24 24 

and three natural pairs of linear equations: 

i/o, ,• iTt + Aq 12m + e 12m~ 

[a = a ' h, h = — — r), (y = — r, r = — — — y), (x — -, e = ; x). (62) 

v ' 3A 2 h vy 3A 2 b+ yh v 6 b- ' v 

The last two pairs, (y,r) and (x,e), are clearly harmonic oscillators. To analyze these 
equations further along with the first pair, we introduce the new variable combination 

z ■= m + y + A h. (63) 

We then find three late time normal modes: 

x + uj 2 _ x = 0, where uj 2 _ = — — (64) 

6 

9 o 4A m + 

y + u+y = 0, where uj + = (65) 

z =0. (66) 

Reexpressed in terms of the late time normal modes the late-time energy constraint is 

, 3 2 (3A 2 , o b + 2 \ / _ 2 6~ 2 \ 

const. = —a Kp = —7-z + I — —m y H — -r I + I 3m x H — -e I , (67) 
Aq \^Aq 24 / \ 24 / 

where each bracket is constant. The physical and geometric significance of two of the normal 
modes is clear, since they directly correspond to the two torsion magnitudes (alternately the 
associated curvature scalars). Thus for this model at late time we find that the two dynamical 
connection modes are essentially dynamical torsion modes and the description "torsion 
cosmology" is phenomenologically appropriate, although it could lead to a misapprehension 
as to the true fundamental dynamical fields @ 

The remaining mode z corresponds to a certain combination of the frame/metric scale 
expansion factor and the + torsion: 

Z = m + f + A H (68) 

which evolves according to 

Z = m- X 2 - m + f(f + 3H) - 2A H 2 + ^. (69) 

o 

§ Although the terminology "torsion cosmology" has often been used in the past to describe the sort of model 
we are considering here, it has recently come to our attention that such terms are fundamentally inappropriate 
and can inhibit a deeper understanding. The basic dynamic variables in this theory are the orthonormal frame 
and connection. 
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(With a nonvanishing cosmological constant this equation would pick up an extra 2 A/3 term.) 
From the above we see that at late time (with vanishing cosmological constant) the Hubble 
expansion rate has the form 

H = a~ 3/2 x const. - ^-/, (70) 

with the + torsion amplitude / oscillating at the frequency uj + . Only the + mode effects the 
expansion rate at late times. The late time acceleration is 

a = a~^^-r, (71) 

which has periodic oscillations at the rate cu + . In this model sometimes the expansion rate is 
accelerating and sometimes it is slowing down. 



5.1. Numerical test 

The validity of our late time analytic results has been tested numerically. Taking the 
parameters as A — 1, A 2 — 0.23, A3 = —0.35, b + = 1.1, and b~ = 0.3, we find uj + = 4.20 
and co>_ = 1.4. They have the relation uj + = 3u;_. Using the same parameters, we plot a full 
and a linear asymptotic normal mode evolution of all the 6 dynamical equations. The behavior 
of the normal mode equations has been observed with several sets of initial values. 

We have plotted one typical case in Fig. [IJ Here we show at late time the asymptotic 
amplitudes: first, the Hubble function, h, second, the normal mode combination of the torsion 
and Hubble function, z, third, the spin-0 normal mode, x, fourth, the spin-0 + normal mode, 
y. The (black) dashed lines represent the exact evolution and the (red) solid lines represent 
the late time asymptotic normal mode behavior. 

As expected, we found that the late time equations are indeed a good approximation. 
The z(t) function approaches a flat line at late time. The plots show that the frequency and 
amplitude of the full and the linear approximation + and 0~ asymptotic equations are very 
close, although there are apparently still some nonlinear effects. 



6. Numerical Demonstration 



In this section, we present the results of a numerical evolution of our cosmological model. 
For all these calculations we take A = 0. Since there is one scalar mode and one pseudoscalar 
mode in this model, it is natural to investigate the interaction between these two modes. We 
find that the pseudoscalar connection mode can generate the scalar connection mode, but not 
conversely. This offers a reason for the existence of the scalar mode, since it is believed that 
the pseudoscalar mode exists and plays an important role in the early universe due to its direct 
interaction with matter. In Section I6.2L we first extend the numerical demonstration of the 
earlier work [2] by including the spin-0~ mode, and then compare our numerical results with 
the observational supernovae data. Not surprisingly we find that the supernovae data can be 
better fitted with the two- scalar-connection-mode model. 
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Figure 1. Asymptotic evolution of the Hubble function, h, the late time normal mode 
combination of torsion and Hubble function, z (note the scale), the late time spin-0~ 
normal mode, x, and the late time spin-0 + normal mode, y. The (black) dashed lines 
represent the actual late time asymptotic evolution and the (red) solid lines represent the linear 
approximation normal modes. 



We need to look into the scaling features of this model before we can obtain the sort of 
evolution results we seek on a cosmological scale. In terms of fundamental units we can scale 
the variables and the parameters as 



^3 -> A 3 , 



E 
A 



2 A, 



(72) 



t — > t/£, a — > a, H 

A ^A , A 2 ^A 2 , A 3 ^A 3 , b + ^b + /f, b- 

where £ 2 = k = 8tcG. So the variables and the scaled parameters b + and b~ become 
dimensionless (from the Newtonian limit A = 1). Equations. (1241 - 1291 ) remain unchanged 
under such a scaling. However, as we are interested in the cosmological scale, it is practical 
to use another scaling — mathematically to make the numerical values of the scaled variables 
less stiff for the numerical integration, and physically to see changes on the scale of the age 
of our Universe. In order to achieve this goal, let us introduce a dimensionless constant T , 
which represents the magnitude of the Hubble time (T = H Q l = 4.41504 x 10 17 seconds). 
Then the scaling is 

t-Tnf, a^a, H — > H/T , f^f/T , X - x/T , R - R/T*, E — > E/Tq, 



oh 
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A -^A , A 2 ^A 2 , A 3 ^A 3 , 6+^T 2 6+, r^T V, A - A/T 2 . (73) 
With this scaling, all the field equations are kept unchanged while the period T — > T T. 

6.1. The interaction between the scalar and pseudoscalar mode 




\ i 
\ i 
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Figure 2. Evolution of the components of the torsion. The top panels show the evolution 
for nonvanishing / and R, and vanishing \ and E initially, corresponding to Case (a). The 
bottom panels show the evolution for nonvanishing \ an d E, and vanishing / and R initially, 
corresponding to Case (b). 

To understand the interaction of these two propagating scalar connection modes, we 
consider two different situations: (a) the pseudoscalar mode vanishes initially with a 
non-vanishing scalar mode; (b) the scalar mode vanishes initially with a non-vanishing 
pseudoscalar mode. The parameters and the initial values for Case (a) are set as 

A 2 = 0.5, A 3 = 1, b + = 2, b~ = 1, (74) 

and 

/(0) = -0.3, #(0) = 0.3, x(0) = 0, E(0) = 0, (75) 

where a(0) = 50 and H(0) = 1 in both the cases. The parameters and the initial values for 
Case (b) were chosen as 

A 2 = 1.0, A 3 = -0.1, b + = 1.5, b~ = 1, (76) 



Cosmological dynamics with propagating Lorentz connection modes of spin zero 



17 



and 

/(0) = 0, R(0)=0, x(0) = 0.3, £(0) = 0.3. (77) 

The numerically calculated evolution of these two cases are shown in Fig. |2l Case (a), which 
corresponds to the top two panels of Fig. [2l shows that the scalar mode cannot generate an 
initially vanishing pseudoscalar mode. It is clear that x an d E stay zero with a dynamic 
spin-0 + mode. The two bottom panels of Fig. [2l which correspond to case (b), show that 
the pseudoscalar mode can generate an initially vanishing scalar mode. It is known that the 
pseudoscalar connection mode will couple to elementary spinning particles. One can expect 
that this mode will be excited by spinning particles in the early Universe, since there could be 
a high spin density during this epoch. Once the pseudoscalar connection mode is generated, 
the scalar connection mode will be excited through the interaction with the pseudoscalar 
connection mode. 



6.2. Accelerating universes 

We would like to compare the numerical evolution values for this model (with A = 0) with 
the observational data of our Universe. The Hubble constant at present, H(t = 1), is 

H = 1 I ~ 70 km (m 

4.41504 x 10 17 s ~ s-Mpc' 

The initial data is set at the current time to = 1, and current value of the Hubble function is 

scaled to unity in this work, just as in [0. The parameters and initial conditions chosen for 

our first case are as follows: 

A 2 = 0M, A 3 = -0.35, b + = 1.1, b- = 0.091, (79) 

and 

a{to = 1) = 50, H(to = 1) = 1, /(to = 1) = -0.335, 

x (t = 1) = 0.378, R(t = 1) = 2.18, E(t = 1) = 2.21. (80) 

The results of the evolution with these parameters and initial conditions are plotted in Fig. [3] 
The expansion factor a is plotted in the top-left panel. In the top-right panel the Hubble 
function H is damped-oscillating at late time. In the bottom-right panel, it is obvious that a is 
damped and oscillating during the evolution and is positive at the current time which 
means the expansion of the universe is currently accelerating. The torsion and curvature 
scalars, f(t), x{t)> R{t), and E(t), are also plotted in the middle panels of Fig.|3]to show the 
correlation of the evolution between these variables. We observe that the frequencies of the 
pairs (x,E) and (f,R) are usually different. The behavior is consistent with the analysis in 
Sec. |5J In order to have a deeper understanding of the settings of this case, the matter density 
p and the effective mass density of the dynamical connection p r are plotted in the bottom-left 
panel. The value of p is decreasing as the universe is expanding and is always positive with 
pa 3 =const, while pr, plotted in the same panel, shows a "damped-oscillating" behavior. The 
damped-oscillating behavior of pr simply indicates that the effective energy density p r is not 
positive-definite in general. We also plot the phase diagrams for Case I in Fig. @]to show that 
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Figure 3. Evolution of the expansion factor a, the Hubble function, H, the scalar and the 
pseudoscalar torsion components, / and x, the affine scalar curvature, R, the pseudoscalar 
curvature, E, the mass densities, p and pr, and the 2nd time derivative of the expansion factor, 
a, as functions of time with the parameter choice and the initial data in Case I. 
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Figure 4. The phase diagrams for Case I. The phase diagram of (F, \, a) is shown in the left 
panel. The (red) solid line is the trajectory of the (F, x, a) evolution starting from the initial 
value (—0.335, 0.378, 50). The (gray) dashed line is the convergence line (0, 0, a) for this 
diagram. The phase diagrams of (F, H, R) and of (x, H, E) are shown in the right panel. 
The (red) solid line is the trajectory of the (F, H, R) evolution starting from the initial value 
(—0.335, 1, 2.18), the (blue) dashed line is the trajectory of the (\, H, E) evolution starting 
from the initial value (0.378, 1, 2.21), and the (black) filled point marks the asymptotic focus 
point (0, 0, 0). 

Table 1. The initial data and parameters for cases I, II, and III. Here the parameter Aq = 1, 
a(t = 1) = 50, and H(t = 1) = 1 in all three cases; t = 1 means t = now, under the scaling 
ofEqs. d72H73]l. 



Case 


A 2 


^3 


b + 


b~ 


/(I) 


x(i) 


R(l) 


E(l) 




I 


0.83 


-0.35 


1.1 


0.091 


-0.335 


0.378 


2.18 


2.21 


0.23 


II 


0.52 


0.475 


1.05 


0.35 


-0.318 


0.225 


2.7 


-1.2 


0.29 


III 


0.635 


0.5 


1.06 


0.42 


-0.361 


0.058 


2.442 


-1.8 


0.27 



the orbit of (f,x> a ) is convergent to a line (0,0,a), and the orbits of (f,H,R) and of (x,H,E) 
both converge to the point (0,0,0). 

In this case the scaled value of p[t = 1) = 0.68 and its physical value is pit = 
T ) = 2.15 x 10~ 30 g/cm 3 . The Universe is supposed to be very close to the critical density, 
p c = 3c 2 H 2 /8nG = 9.47 x 10~ 30 g/cm 3 ; we find the ratio Q m = p/p c w 23%. In the standard 
ACDM model, f2 m ~ 30% with 5% baryonic matter and 25% dark matter. For our model 
Vl r = pv/ p c = 77% acts like the energy density of dark energy. Therefore, this dynamic 
connection model is able to describe a presently accelerating expansion of the Universe with 
a proper amount of matter density. From the field equations we can see that the effect of the 
"dark energy" mainly comes from the nonlinearity of the field equation driven by the dynamic 
scalar and pseudoscalar connection modes. We also found other cases, two of which are listed 
in Table CD along with Case I; they are obtained by taking different values for the parameters 
and the initial conditions. We find that the results for the other two cases have a behavior 
qualitatively similar to that of Case I. 



Cosmological dynamics with propagating Lorentz connection modes of spin zero 20 




Figure 5. Comparison of different spin-zero connection models and the standard ACDM 
model with the observational data via the relation between the distance modulus /i and the 
redshift z. The supernovae data points, plotted with (brown) circles, come from [ 38 1 . The 
result of the standard ACDM model (O m = 0.3, = 0.7) is plotted by the bold solid 
line. The results of Case I, II, and III are represented by the (red) dashed line, the (green) 
dot-dashed line, and the (blue) dotted line, respectively. In the inset, the models and data are 
shown relative to an empty universe model (£1 = 0). 



We compare our results with the supernovae data. Distance estimates from SN la light 
curves are derived from the luminosity distance 

where L- mt and T are the intrinsic luminosity and observed flux of the SN, and the redshift 
z = a{\)/a{t) — 1. Logarithmic measures of the flux (apparent magnitude, m) and luminosity 
(absolute magnitude, M) were used to derive the predicted distance modulus 

/i = m-M = 51og 10 d L + 25, (82) 

where m is the flux (apparent magnitude), M is the luminosity (absolute magnitude), and 
(II in the formula should be in units of megaparsecs. We found the relations between the 
predicted distance modulus /i and the redshift z in the three cases; they are plotted in Fig. [5] 
For comparison, we also plot the prediction of the ACDM model with f2 m = 0.3 and = 0.7 
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by employing the following formula Il38l 



y/(l + zf{l + Q m z) - z{2 + z)Q A ' 
The astronomical observational data 113811391 are also plotted in Fig. [5] for comparison. The 
plots show that for small redshift z (e.g., z < 1.9) all three cases of the dynamical connection 
models give an accelerating universe just like the ACDM model does. For larger z these cases 
might turn the Universe into a deceleration mode, which is consistent with the behavior of the 
various quantities shown in Fig. [3] We can see that Case I gives the closest curve behavior 
to the one from the ACDM model. In Fig. |5l we demonstrate the possibility of the spin-zero 
connection fields accounting for the effect of dark energy with a suitable set of parameters and 
initial data. A comparison of Fig. [5] with the results in [2] shows that this two-scalar-mode 
model can give (not surprisingly) a better fitting of the supernova data than the one-scalar- 
mode model can. 



cT (l + z) 



7. Discussion 

From a series of earlier works lfT6l IT5l [P71 [T8l it was concluded that the Poincare Gauge 
Theory of gravity has two good dynamic Lorentz connection modes, the "scalar" mode (spin 
+ ) and the "pseudoscalar" mode (spin 0~) which satisfy 2nd order equations. 

Here we extended a previous work 0, which considered a PGT cosmological model 
with one dynamic Lorentz connection mode having spin + , to the case where both the 
scalar and pseudoscalar connection modes are dynamic. The objectives are (i) to study 
this PGT cosmological model (and in particular how well it can match the present universe 
observations) and (ii) to get a deeper understanding of the dynamics of the PGT. 

From the cosmological homogeneous and isotropic assumptions the scalar and 
pseudoscalar curvatures R, E and the temporal components of the trace torsion / and axial 
torsion x survive and affect the evolution of the universe in this two-connection-mode model. 
Recognizing the equivalence of the model to one describing a particle with three degrees 
of freedom, we constructed an effective Lagrangian and the corresponding Hamiltonian by 
imposing the FLRW symmetry on the field theory action. The system of ODEs obtained 
therefrom are the same as the evolution equations obtained by imposing the FLRW symmetry 
on the equations derived from the PGT Lagrangian density. 

With the evolution equations (l24)) - (|29l) and the associated energy constraint (|3Q|) we 
analyzed the late time asymptotic expansion. We found there are three normal modes: one 
related to the Hubble expansion and two dynamic modes represented by the torsion and 
curvature components. It was found that only the scalar mode affects the late-time expansion 
rate. The numerical analysis focused on the interaction between these two modes, the study 
of the possible behavior of the Universe, and the fitting to the observed supernova data. It 
was shown that the dynamical activity of the pseudoscalar mode could excite the scalar mode 
via the nonlinear coupling of these two modes, but the converse does not happen: one can 
have the scalar mode excited without any pseudoscalar excitation. Like the one-mode model 
in 0, the present model allows for an expanding universe with an oscillating component in 
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the expansion rate. Consequently, although on the average the expansion is slowing down, 
the universe can have an accelerating expansion at the present time. The additional degree 
of freedom in this two-mode model compared to the one in J2||, not surprisingly, allows us to 
obtain a better fit to the supernova data. 

From the evolution equations (|26l) -(T29l). we can see the nonlinear coupling between the 
scalar and pseudoscalar modes. In Sec. 16.11 we demonstrated the excitation of the scalar 
mode through the dynamical activity of the pseudoscalar mode. Such a nonlinear interaction 
between these two modes offers a natural mechanism to fuel the strength of the scalar mode 
in the evolution of the Universe. We stress that there is no known fundamental material 
source which directly excites the scalar mode; this part of the Lorentz connection simply does 
not interact in any obvious fashion with any familiar type of matter [12111 . Conversely, the 
pseudoscalar mode is naturally driven by the intrinsic spin of fundamental fermions; in turn it 
naturally interacts with such sources. Indirectly, the + mode could be enhanced and activated 
dynamically through the aforementioned non-linear mechanism, in addition to any possible 
primordial amplitude from the early universe. 

From the late-time analysis in Sec. [51 we showed that only the scalar mode, not the 
pseudoscalar mode, plays a direct role in affecting the expansion rate of the Universe. This 
result is perfectly consistent with our understanding of the characteristics of these two modes: 
Due to the ability of interacting with fermionic matter, it is generally thought that the axial 
torsion (controlled by the pseudoscalar part of the Lorentz connection) must be small and 
have small effects at the present time |fT9l . Conversely, the scalar Lorentz connection mode 
could be considered as a "phantom" field, at least in the matter-dominated epoch, since it will 
not interact directly with matter, and yet can drive the Universe in an oscillating fashion with 
an accelerating expansion at the present time. 

As discussed in [2], the two Lorentz connection spin-zero modes in this model are in 
some ways effectively like a scalar field and a pseudoscalar field, yet these two "scalar" fields 
are fundamentally different from the various scalar field models of unknown matter, e.g., the 
quintessence models, in the following ways: 

• this cosmological model is derived naturally from a geometric gravitational theory, the 
PGT, which is based on fundamental gauge principles, instead of on the hypothesis of 
the existence of a dark energy tailored to producing an explanation of an accelerating 
universe; 

• there are, consequently, only a few free parameters in this cosmological model, instead of 
an ad hoc potential that can be rather arbitrarily chosen to fit the observations. Therefore, 
this PGT cosmological model should be more restrictive, and should be easier to be 
confirmed or falsified; 

• based on its geometric character, the coupling of the dynamic parts of the Lorentz 
connection to the other fields is nothing like that which has ever been advocated for 
hypothetical scalar fields. 

Thus this PGT cosmology with a Lorentz connection having dynamic "scalar" modes and the 
quintessence models are characteristically different, even though there are some similarities. 



Cosmological dynamics with propagating Lorentz connection modes of spin zero 



23 



As mentioned in our previous work, if we consider the spacetimes as Riemannian instead 
of Riemann-Cartan, by absorbing the contribution of the post-Riemannian terms of this model 
into the stress-energy tensor on the rhs of the Einstein equations, as indicated in (I33I34I35I36I) . 
then this contribution will act as a source of the Riemannian metric, effectively like an exotic 
fluid with its mass density pr and pressure pr varying with time (although the time evolution 
of these torsion and curvature terms are not like that of any fluid). Moreover, the effective 
fluid will appear to have presently a negative pressure, and consequently a negative parameter 
in the effective equation of state, i.e., ur = pv/ Pr, which drives the universe into accelerating 
expansion. Note that there is no constraint on the value of cur which appears here, and its 
value could vary from time to time. It should be stressed that this is not a real physical fluid 
situation; the truth is that c^r is nothing like "a connection field equation of state", it is just 
a proportionality factor between pr and pr, two expressions which effectively summarize the 
contribution of the connection (via the curvature and torsion) acting as a source of the metric. 
The ratio cu r is of interest only to help understand the acceleration of this model and to enable 
a limited comparison with other dark energy proposals. 

By imposing the FLRW symmetry on the Lagrangian density, we constructed an effective 
Lagrangian as well as the corresponding Hamiltonian. One benefit of the former is a simpler 
derivation of the dynamic equations (|24|)-(|29|). The latter should also prove useful, as 
the Hamiltonian formulation is the framework for the most powerful known techniques for 
analytically investigating the dynamics of a system. By these techniques one can to apply 
the experience accumulated in dealing with conservative classical mechanical systems. An 
effective Lagrangian and the corresponding Hamiltonian allows one to visualize the system 
as a particle moving in a potential. This would be very helpful in gaining a better appreciation 
of the dynamics of any sophisticated model. (Note, it is not necessarily true that an effective 
Lagrangian can be found in an arbitrary cosmological model. Extrapolating from GR, one 
can conjecture that this is possible for all PGT Class A Bianchi models with suitable sources, 
including pressure and spin). As we have seen in Sec. [51 the effective mechanical system 
methods were also useful for the late time normal mode analysis. 

There have been some studies on PGT cosmology with dynamic scalar connection modes 
since the model proposed in tTJ|2]|. Wang and Wu 113611 considered a related, but fundamentally 
different model, which turns out to have only the dynamic CT mode. They considered the early 
universe and showed how in their model such a dynamic PGT connection could account for 
inflation (for another approach to using the PGT to account for inflation, see OTI ). Li et 
al. 01 HI presented a nice analysis of the scalar mode model of from a more mathematical 
angle in order to get a deeper insight into the behavior of the dynamical system. In their work, 
they found the critical points of the system and the corresponding ranges of the parameters. 
In the latter work they also fit the model to the supernova data to find the best fit values of 
the parameters. These works considered quite general ranges of the parameters and found 
several interesting dynamical effects. We note that many of these interesting effects happen 
in parameter ranges which are outside of the restrictions considered by [2J to be physically 
necessary in order to have good linear modes (long ago (H [13J the conditions were found so 
that the propagating modes should carry positive energy and satisfy the no-faster-than-light 
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condition). Also, a matter density of about 25% of the critical density was also imposed in 
to give a more physical meaning to the a curve fitting. Further investigations and a careful 
study of the model will be needed to pin down the acceptable ranges of the parameters. In the 
present work, we have chosen the range of the related parameters following the result of ifTTl 
for good propagating scalar linear modes, just as in 0. Applying the methods used in ifJUH 
to the present extended model would surely lead to further insights. 

One may wonder: how large must the post-Riemannian fields be in order to produce 
observable effects in the the present day universe, e.g., the observed acceleration? Conversely, 
how large can the torsion or curvature scalars be without violating some observational 
constraint? The questions merit a detailed study. Here is a simple argument that indicates 
a magnitude. Let us compare the terms in the Lagrangian density and the field equations for 
the model in which the PGT Lorentz connection has scalar dynamical modes and the Einstein 
theory with a cosmological constant. (In our present work we have deliberately included in 
most of the dynamical equations a possible cosmological constant; this was done not only for 
greater generality but also to facilitate just such a comparison. In our numerical evolution for 
our model we used A = 0.) Note that the presumed cosmological constant is "so small" that 
it has no noticeable effect in the laboratory, nor on the solar system scale, nor on the galactic 
scale. Nevertheless it is large enough to have the dominant effect on the cosmological scale. 
Hence we are led to infer that we should consider that one or more of the post-Riemannian 
terms (A 2 f 2 , A 3 x 2 , b + R 2 , b~E 2 ) should be comparable to the cosmological constant (which 
is about 3p) in the ACDM model. With such a choice we can expect that the post-Riemannian 
terms may be able to accelerate the universe and yet not be conspicuous on smaller scales. 

The introduction of a new ingredient (i.e., the CT connection mode which is reflected 
in the axial torsion and the pseudoscalar curvature) in this work raises the concern of 
the experimental and observational constraints on this field. There have also been some 
laboratory tests in search of torsion Il40l |4~U| . The main idea among these experiments is 
the spin interaction between matter and torsion. The theoretical analyses and the high energy 
experimental data on four-fermion vertices sets the lower bound for the (pseudoscalar) torsion 
mass > 200 Gev [fT9l l20l I2T1 142| . The cosmological tests on torsion have investigated 
the effect of torsion-induced spin flips of neutrinos in the early Universe — which could 
alter the helium abundance and have other effects on the early nucleosynthesis [|43l |44"1| . 
From Table [Q the parameters chosen for the range of the the torsion mass are consistent 
with the aforementioned analyses. Our model is also comfortable with the most restrictive 
experimental limits found on torsion [|45l . For torsion being applied to the cosmological 
problem, Capozziello et al. Il35l l46l have done a serious study on replacing the role of the 
cosmological constant in the accelerating Universe. With a totally antisymmetric torsion 
without dynamical evolution, their model is consistent with the observational data by tuning 
the amount of the torsion density. Compared with them, the model in this work allows the 
pseudoscalar torsion to evolve dynamically. This difference might enable a more fruitful 
physics to be studied. 
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8. Conclusion 

In this work we considered the two "scalar" dynamical modes of the PGT Lorentz connection 
in a cosmological setting and have proposed it as a viable model for explaining the current 
status of the Universe. Besides seeking a better understanding of the PGT, we have considered 
the prospects of accounting for the outstanding present day mystery — the accelerating 
universe — in terms of an alternative gravity theory, more particularly in terms of the PGT 
with a dynamic Lorentz connection having only two dynamic modes, carrying spin-0 with 
even and odd parity. With the usual assumptions of isotropy and homogeneity in cosmology, 
we find that, under the model, the Universe will have with generic choices of the parameters 
an expansion rate which oscillates. The connection in this model could play the role of dark 
energy. With a certain range of parameter choices, it can account for the current status of the 
Universe, i.e., an accelerating expanding universe with a value of the Hubble constant which 
is approximately the present one. Thus we have considered the possibility that a certain 
geometric field, a dynamic Lorentz connection — which is naturally expected from spacetime 
gauge theory — could fully account for the accelerated universe. 

The + mode, which directly drives the acceleration of the universe, does not couple 
directly to any known material source. By way of non-linear terms it could come indirectly 
from the huge density of the particles with sufficient spin alignment in the early universe 
which directly excite the 0~ connection mode. The + mode could be considered as a 
"phantom" field, at least in the matter-dominated epoch, since it will not interact directly 
with matter; it only interacts indirectly via the gravitational equations. Then the dynamics of 
the scalar torsion mode could drive the Universe in an oscillating fashion with an accelerating 
expansion at present. It is quite remarkable that a gauge theory of dynamic geometry naturally 
presents us with such a "phantom" field. This natural geometric field could act like a dark 
energy. 
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Appendix A: The choice of parameters 

Regarding our choice of parameters. From the table which can be found in any one of 
lfl5l[P71[T8ll . we find that to kill the dynamics of the 1 + , 2 + , l~, 2~ modes we want to take, 
respectively, 

b 2 + h = 0, 6i + h = 0, h + h = 0, b 1 + b 2 = 0. (84) 
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For dynamic + and 0~ we want, respectively, 

h + b 6 > 0, and b 2 + b 3 < 0, (85) 

from Table 3 in lfT8l . Now, due to the Bach-Lanczos identity, we can choose any one of the 
parameters b k to vanish. Taking, say b 4 = we then get that we also want b\ = b 2 = b 5 = 0, 
leaving 6 6 > 0, 5 3 < 0. We also find for the dynamic + and 0~ modes, from Table 3 in lfT8ll 
the respective restrictions 

a a 2 (2a + a 2 ) < 0, and a + 2a 3 < 0. (86) 

In terms of the parameters used in the present work, i.e., 

b + = b 6 , b~ = -b 3 , A k = -a k , m + = A + A 2 /2, m~ = A + 2A 3 , (87) 

these restrictions become 

b + > 0, b~ > 0, A A 2 m + > 0, m~ > 0. (88) 

The Newtonian limit gives A — 1. A positive kinetic term in the action requires A 2 > 0. 

In order to facilitate a comparison of the works of various groups, we here include the 
parameter conversion between those of the Cologne group of Hehl and coworkers (which we 
follow), Minkevich and coworkers flM HUSHED EH, and Goenner and Miiller-Hoissen [fTTTl . 
In Hehl's work, the parameters related to the model described in this paper are a , a 2 , a 3 , b 3 , 
b$. Goenner and Miiller-Hoissen used c\ ■ ■ ■ eg as the parameters in their work. The Goenner 
and Miiller-Hoissen parameters are related to Hehl's by 

a 3 a 2 ao 

Cl + 3c 2 = —, a = c 1 + 3c 3 = —, c 4 = -y, 

c 9 - c 8 = -, c = 2(6c 5 + 2c 6 + 2c 7 - c 8 - c 9 ) = — , (89) 



and thus 



, m 
2c4 — a = m , C4 — 2ci — 6c 2 = ——. (90) 



The parameters used by Minkevich are related to Hehl's by 



and thus 



a a b 3 h-b 3 b 6 

as, = a 2, Jo = -y, Q2 = j, qi = — ^ — , / = y , (91) 



2/ + ^ = m + , / -6=— . (92) 
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